Role of the Pauli principle in collective-model coupled-channels calculations 
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A multi-channel algebraic scattering theory, to find solutions of coupled-channel scattering prob- 
lems with interactions determined by collective models, has been structured to ensure that the Pauli 
principle is not violated. By tracking the results in the zero coupling limit, a correct interpretation 
of the sub-threshold and resonant spectra of the compound system can be made. As an example, 
the neutron-^^C system is studied defining properties of ^''C to 10 MeV excitation. Accounting for 
the Pauli principle in collective coupled-channels models is crucial to the outcome. 
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At energies above 25 MeV, by using optical potentials 
formed by full folding effective two-nucleon interactions 
with microscopic (nucleon based) descriptions of the tar- 
get structure, the importance of treating the Pauli princi- 
ple has been well established However, in the domain 
of low-energy nucleon scattering for which an explicit 
coupled-channels theory of scattering is essential, the sig- 
nificance of Pauli exclusion effects has not been well de- 
fined. Many coupled-channels codes are available, some 
of which perform phenomenological collective-model cal- 
culations searching on parameter values of the chosen 
function forms to find a best fit to experimental data. 
But while it has long been known that any such models 
violate the Pauli principle 0, 1^, quantification of that 
violation is lacking. 

To study the effects of the Pauli principle in a macro- 
scopic (collective model) approach is not a trivial task. 
In a recent publication M, the orthogonalizing pseudo- 
potential (OPP) method Q was generalized to treat 
this problem. That was a small though important part of 
the full theoretical framework of the multi-channel alge- 
braic scattering (MCAS) theory of scattering 4J. Therein 
the OPP was used in finding the spectra, bound and res- 
onance properties, of ^'^C. However, implications of the 
role of the Pauli principle in collective model coupled- 
channel calculations arising from the use of the OPP 
was not discussed. Such is a purpose of this letter. 
Another is that the method could be pertinent for any 
study requiring coupled channel solutions of quantal sys- 
tems involving fermions. As the example, we study the 
effects introduced by the Pauli principle in collective, 
geometrical-type, models for low-energy nucleon-nucleus 
processes that can be characterized from the spectrum 



of the compound nucleus. That spectrum includes the 
states that lie below the nucleon-nucleus threshold and 
in the continuum as revealed by the narrow and broad 
resonances that lie upon a smooth but energy dependent 
background of the elastic scattering cross section. This 
can be done in a systematic and self-consistent way since 
the MCAS approach facilitates such a determination of 
the sub-threshold bound states and resonances of the 
compound nucleus. This theory, with which one solves 
the coupled-channel Lippmann-Schwinger (LS) equations 
for the nucleon-nucleus system considered, is built upon 
sturmian expansions of an interaction matrix of potential 
functions. 

The MCAS method has been developed in momen- 
tum space and the starting matrix of potentials may 
be formed by folding effective two-nucleon interactions 
with one-body density matrices of the target (studies in 
progress) or, as is more common, from a collective model 
description of the target states and excitations. As in 
that recent publication 3i have used a rotational 
collective-model representation with deformation taken 
to second order. We chose Woods-Saxon functions and 
their various derivatives to be the form factors for all 
components each with characteristic operators of diverse 
type. The interactions were allowed to depend on parity 
as well. With such a characterization, we were able to 
describe all important aspects, at positive and negative 
energies, in the neutron-^^C system. 

With the MCAS approach and a collective model pre- 
scription for the starting matrix of potentials, the OPP is 
used in the process by which the sturmians are specified. 
The OPP inclusion ensures that all sturmians in the (fi- 
nite) set selected as the basis of expansion of the matrix 
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of potentials contain few or no components equivalent to 
the external nucleon being placed in an already densely 
occupied orbit. That scheme is an approximation as we 
discuss later by assessing the spectra of 12,i3q g^j^j ^j^g 
single neutron spectroscopic amplitudes that link them 
using information obtained from large space no-core shell 
model calculations. But it is a good approximation. 

The role of the Pauli principle is studied by comparing 
results found with and without using the OPP scheme to 
select the sturmians that form the expansion set. Note 
that the actual matrix of potentials is the same through- 
out though extra information on single nucleon plus a 
core nucleus state underlying each sub-threshold bound 
and resonance in the compound system has been ob- 
tained by taking the zero deformation limit. 

Full details of the MCAS scheme have been published 
and the reader is referred there for those, as well as 
for specifics of the notation we use herein. In momentum 
space for potential matrices Vcd {p, q) , one seeks the so- 
lution of coupled LS equations [see Eq. (1) in Ref. 0|], 
which involve both open and closed channel contribu- 
tions. With incident energy E, the channel wave numbers 
for the open and closed channels are kc and he respec- 
tively. Solutions of those LS equations are sought using 
expansions of the potential matrix elements in (finite) 
sums of energy-independent separable terms, 
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Vcc'(p,q) ~ ^Xcn{p)Vn^Xc'n{q) , 



(1) 



where Xc'n{q) are the Fourier-Bessel transforms of the 
selected sturmians whose eigenvalues are To predict 
observables one requires the multichannel S'-matrix. In 
terms of the multi-channel T-niatrix, that has closed al- 
gebraic form 
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Tec - Xc„(fcc)([r7-Go]-i)„„, xc'«'(fcc') ,(2) 

n,n' — 1 

where now c, c' refer to open channels only. In this rep- 
resentation, Go and f] have matrix elements 
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The bound states of the compound system are defined by 
the zeros of the matrix determinant when the energy E 
is negative, and so link to the zeros of {\r] — Gq\} when 
all channels in Eq. are closed. 

As noted above the sturmians are solutions of homoge- 
neous Schrodinger equations for the matrix of potentials. 



TABLE I: Shell occupancies of protons (or neutrons) in states 
of i^c. 



orbit 


0^ 


2^ 




Osi 

2 


1.963 


1.962 


1.968 


Op 3 

2 


3.054 


2.858 


3.075 


Opi 

2 


0.842 


1.028 


0.804 


higher orbits 


0.124 


0.129 


0.120 



In coordinate space if those potentials are designated by 
local forms Vcc'{r)S{r — r'), the OPP method is to use 
sturmians that are solutions for nonlocal potentials 



Vac'ir, r') = Vac'{r)Sir - r') + \A,{r)A,{r')Sc 



(4) 



where A{r) is the radial part of the single particle bound 
state wave function in channel c spanning the phase space 
excluded by the Pauli principle. The OPP method holds 
in the limit A — > oo, but use of A = 100 MeV suffices. 

The spectrum of ^^C also was calculated in the shell 
model using the program OXBASH and with the 
MK3W interaction Q. The positive parity states of -^^C 
were calculated in a complete (0 -I- 2)Huj space using this 
interaction, while the negative parity states were calcu- 
lated in a restricted {l + 3)huj space. In both calculations 
the same single particle basis of Os up to and including 
the 0/lp shell was used. Hence the restriction from a 
full (1 -I- 3)huj study is that we have not included the 
0gld2s shell. With exceptions, most notably the super- 
deformed state at 7.654 MeV and the known collec- 
tive 3^ state at 9.64 MeV, the calculated spectrum to 
20 MeV excitation agrees well with observation 0. So 
also do results of calculations of elastic and inelas- 
tic scattering data (form factors from electron scattering 
and differential cross sections and analyzing powers from 
proton scattering) without the need for any a posteriori 
core polarization corrections. 

Of interest here are the details of the low lying spec- 
trum. First, in Table |I| we list the nucleon shell occu- 
pancies in the three lowest states of ^^C. Clearly the Os- 
and Op- shells have dense occupancy: essentially 4 nucle- 
ons filling the Os shell while there are almost 8 nucleons 
in the Op shell. Those eight nucleons are distributed be- 
tween the sub-shells, so blocking the Opa orbit as we do 
in using the OPP method is an approximation. Note also 
that the second excited state of ^^C, the 0^(7.654 MeV) 
is well known to be a super-deformed 3a chain and, as 
such, a much larger space is needed for a good descrip- 
tion. The lowest three 0+ states in our shell model are 
the ground, at 12.25, and at 23.03 MeV. They are much 
more spread than measured energies and have structure 
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TABLE II: Dominant components of shell model spectro- 
scopic amplitudes. Energies in brackets are in MeV. 



TABLE III: Parameter values of the base potential (in MeV). 
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12 c 
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(g-s.) 
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-0.7285 


Op 3 
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-1.0040 


Opi 
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-0.4738 


2 + 
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(3.09) 


Isi 
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-0.9088 


Ods 
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-0.3162 


Isi 
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-0.0605 


3- 
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yo.Oo ) 


n-n ^ 
\jy 3 
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U.40U4 


Oz? 3 
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1 nn/in 
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0/| 


-0.8342 
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(3.85) 


Ods 
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0.8129 


Off 5 
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0.4799 


Ods 
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0.0096 
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0.0840 
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(6.86) 
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-0.2147 


2 


0.5372 


Ods 
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-0.0102 
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-0.0907 
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-0.7714 






5 + 


(8.88) 




-0.0349 


Off 5 


-0.2568 


Ods 


0.2829 


2 
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0^3 
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-0.2694 
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-0.2391 














2 











K)(7r) 


V«<?(7r) 




Vls(jl) 


TT = — 


-49.144 


4.559 


7.384 


-4.770 


TT = + 


-47.563 


0.610 


9.176 


-0.052 


Geometry 


Ro = 3.09 fm 


a = 0.65 fm 


/32 = -0.52 
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Notice that the first dominantly 2fvjj state lies at 
23.03 MeV excitation; a calculated energy that can be 
expected to fall with the addition of higher hw compo- 
nents. That has not been seen sufficiently at least to the 
^hhj level vifith an ah initio shell model So while the 
0;^ state may be the one that is super-deformed, a con- 
vergence in energy will require a greatly increased space. 
But as the 0^ state is not very important in the forma- 
tion of resonances and bound states 0, use of our shell 
model should suffice for the comparisons we make. This 
may have more bearing on future results as we move to 
use a microscopic MCAS in which the matrices of in- 
teraction potentials will be formed using nucleon-based 
structure. 

Next we consider how each state in the low excitation 
spectrum of ^'^C maps onto a single neutron added to any 
of the three selected states of ^^C. The relevant one-body 
spectroscopic amplitudes for /'^ — > J'^ , 
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are listed in Tabled The shell model calculations gave 
more values for addition of that neutron in higher shell 
states, but those spectroscopic amplitudes (not listed) all 
have magnitudes less than 0.1. 

Results of calculations of the neutron- ^^C system re- 
ported previously Q , used the parameter values that are 
specified in Table Hill 

In Fig. n the results are compared with data, both 
elastic scattering cross sections as well as the spectrum 
of ^'^C. Therein it is clear that the three states of ^^C 
suffice to deal with information to ~ 10 MeV excitation 



FIG. 1: The spectra of i^-^-^C and the elastic cross section 
(barn) for n-|-^^C system. The ENDF data (circles) [TO] are 
compared with our full MCAS results (solid line). Note that 
the identification of the ^''C states' spin-parities is 2J'^ . 



in the compound with corroboration in the scattering of 
up to 5 MeV. Spin-parity assignments, bound state en- 
ergies and resonance centroids, widths of the resonances, 
and the background scattering all are very well matched 
by the calculations. A most interesting feature is what 
occurs as the coupling tends to zero. In that limit, all 
of the compound resonances shrink to be bound states 



TABLE IV: Pauli effects on sub-threshold and bound states 
in the continuum in the limit /32 — > (with Vas ~ 0). 





J" 


With Pauli 


No Pauli 


n+i^C 




1 + 
2 




-26.57 


0S1/2+0+ 


3 
2 


+ 5 + 
' 2 




-22.13 


Osi/2 + 2+ 




1 + 
2 




-18.91 


OSl/2 + Ot 




3- 
2 




-8.849 


Op3/2 + 0+ 




1 - 
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-4.685 


-4.685 


Opi/2 + ot 


1-3 
2 ' 2 


-5-7- 
'2 '2 




-4.410 


Op3/2 + 2+ 




3- 
2 




-1.195 


Opi/2 + 0+ 




1 + 
2 


-0.837 


-0.837 


1S1/2+0+ 


3 
2 


- 5 - 
' 2 


-0.246 


-0.246 


Opi/2 + 2+ 




5 + 
2 


-0.171 


-0.171 


0^5/2 + 0+ 




1 - 
2 


2.969 


2.969 


Opi/2 + ot 


3 
2 


+ 5 + 
' 2 


3.601 


3.601 


1S1/2+2+ 


1+ 3 + 
2 ' 2 ' 


5+ 7+ 9+ 
2 '2 '2 


4.267 


4.267 


0^5/2 + 2+ 
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in the continuum. In this hmit, calculations were made 
with the spin-spin interaction strengths set to zero, and 
so offsetting a splitting that is most evident with the two 
odd parity states built from coupling with a Opi neutron. 
The results of these limit calculations are collected in Ta- 
bic IIVI Therein the states are listed in the order from 
most bound to largest continuum energy whether they 
are real or spurious. The energies listed in columns 2 and 
3 respectively were found in the zero deformation limit 
with and without the OPP treatment of Pauli blocking. 
In the last column we display what dominant charac- 
ter (neutron orbit coupled to state in ^^C) is found for 
each state in ^'^C. Disregarding the Pauli principle clearly 
gives many spurious states. However, notice that there 
are matching entries for every resonance state whether 
the Pauli principle is taken into account or not. That 
has led to the erroneous belief that a simple adjustment 
of parameter values is all that is needed to define scat- 
tering cross sections and that the Pauli principle effects 
are unimportant for scattering. Not only is that phe- 
nomenology not guaranteed to work in other cases but 
also the mixing of components caused by finite deforma- 
tion is quite different when the Pauli principle is or is not 
satisfied. A calculation made ignoring the Pauli principle 
gives an incorrect description of all states. 

The resonance centroids tend to three limits. The high- 
est is at 4.267 MeV with five entries from to |^ as 
is formed by attaching a neutron to the 2"*" state 

in ^^C. The second is at 3.601 MeV having two entries 
which equate to a lsi/2 neutron coupled to the 2+ state 
of the target. The third, the only odd parity resonance 
ascertained from these calculations within the range of 
energies to 5 MeV, is identified as a ^ resonance. It lies 
7.65 MeV above the calculated value for the ^'^C ground 
state and can then be associated with binding a Opi/2 
neutron to the second 0"*" state of ^^C. 

The bound states are less clear with regard to domi- 
nant particle coupling character. From shell model cal- 
culations, the ^ (ground state) and the | state are 
sizable mixtures of p-shell nucleon coupling to both the 
ground and 2+ states in '^^C. But the and |^ bound 
states in ^^C are dominantly formed respectively by a 
lsi/2 and a 0^5/2 neutron coupled to the ground state of 
^^C. The energies found in the zero coupling limit sup- 
port the inferences made above. Notably, the bound |^ 
tends to -0.171 MeV in that limit, as 4.267- (-0.171) = 
4.438 MeV; the excitation energy of the 2+ state in ^^C. 
Likewise the doublet | , | tends to 3.601 MeV and as 

the bound state 5^ tends to —0.837 MeV in the zero cou- 
pling limit, 3.601 - (-0.837) = 4.438 MeV, the excitation 
energy of the 2+ state in ^^C once more. Finally the | 
and I states both are bound by —0.246 MeV and so the 
energy gap of that pair from the ^ state (—4.685 MeV), 



is 4.439 MeV. 

But deformation makes significant admixing of these 
nucleon plus core nucleus elements. Indeed it is this spu- 
rious mixing that is the most serious concern about calcu- 
lations that do not take the Pauli principle into account. 
To reveal that, we have repeated the limit calculations 
excising the OPP effects and thereby solving the problem 
in a way equivalent to coordinate-space solutions of such 
coupled channels equations. That calculation clearly has 
given spurious states. 

Summarizing, the MCAS approach has been used to 
evaluate (low-energy) n-^^C elastic scattering and to 
characterize sub-threshold states of ^'^C. A collective 
model prescription with the three lowest states in the 
^^C spectrum was used. The results well match observed 
data but only when allowance for the influence of the 
Pauli principle was made. Without such allowance, many 
spurious states result. Most strikingly, the ground state 
of ^'^C then has the wrong spin-parity and a binding far 
in excess of the known value. But more disturbing is 
that when states may be matched (in energy and spin- 
parity) their underlying nucleon plus ^^C compositions 
are wrong. By tracking results to the /32 — > limit, the 
dominant parentage of each sub-threshold and resonant 
state in this system has been identified. 
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